Intensity Thresholds and the Statistics of the Temporal Occurrence of Solar Flares 
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Introducing thresholds to analyze time series of emission from the Sun enables a new and simple 
definition of solar flare events, and their interoccurrence times. Rescaling time by the rate of events, 
the waiting and quiet time distributions both conform to scaling functions that are independent 
of the intensity threshold over a wide range. The scaling functions are well described by a two 
parameter function, with parameters that depend on the phase of the solar cycle. For flares identified 
according to the current, standard definition, similar behavior is found. 
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The solar corona is a very high Reynolds number tur- 
bulent plasma producing intermittent bursts of radia- 
tion. Plasma forces twist the coronal magnetic fields 
until stresses are suddenly released, an avalanching pro- 
cess governed by magnetic reconnection . The released 
magnetic energy induces radiative emission that can be 
detected as a flare. Flares exhibit scale invariant statis- 
tics. For instance, the probability distribution of flare 
energies is a power law spanning more than eight orders 
of magnitude similar to the Gutenberg-Richter law 

for earthquakes. The distribution of magnetic concentra- 
tion sizes on the photosphere is also scale invariant, and 
the coronal magnetic network embodies a scale-free net- 
work 0,0- In fact, a model of self-organized criticality 
(SOC) with avalanches of reconnecting flux tubes repro- 
duces the observed scale- free network structure 0, |(j . 

As part of the debate on the characterization of mag- 
netohydrodynamic turbulence in this regime 
H, , interest has focused on comparing interoccurrence 
times between flares with those in models of SOC. Analy- 
ses of flare catalogs have indicated scale invariance of the 
waiting times, but the behavior was found to vary with 
the phase of the solar cycle 01 an d with the methods 
used to analyze the catalogs. (See e.g. Ref. 0,0]) The 
prior belief that avalanches occur with Poissonian waiting 
times in the well-known BTW sandpile model ( giv- 
ing an exponential distribution of waiting times) argued 
against the SOC hypothesis ||. However, including a fi- 
nite detection threshold leads to a power law distribution 
of quiet times even for the BTW model , when dura- 
tions and quiet times are measured with the same clock. 
Since the turnover time scale for flux to be regenerated 
in the corona is of the order of ten hours 0], while the 
correlated waiting time intervals between flares can ex- 
tend up to years, the physical mechanism(s) responsible 
for these correlations resides in the turbulent convective 
region beneath the photosphere that generates magnetic 
flux and drives it into the corona. Systematic studies of 



the temporal pattern of flares can give insight into the 
dynamics of magnetic flux in the convective region, or 
the solar dynamo, which is difficult to observe directly. 

Here we show that the interoccurrence times between 
flares has a hierarchical scaling structure when flares are 
defined as intervals during which the emission exceeds 
a given threshold. Rescaling time by the rate of these 
events, we find universal behavior for the interoccurrence 
times, which is independent of threshold. Both at so- 
lar maximum and at solar minimum the scaling function 
can be fitted by a simple two-parameter function. This 
generalized Lorentzian arises naturally within a simpli- 
fied model based on the time-dependent Poisson process. 
From this model one can naturally infer an exponential 
distribution of flaring rates at solar maximum. At solar 
minimum the distribution of quiet- or laminar - times is 
accurately described by on-off intermittency [lfij , a mech- 
anism already proposed to describe the solar cycle [lfij . 

In extremely intermittent time series, like e.g. earth- 
quakes, events are spikes separated by a smooth back- 
ground, and are easily and uniquely defined. This is not 
the case for the solar data analyzed here, where the in- 
tensity decays slowly after a local peak, allowing overlaps 
with subsequent peaks. In this case, the introduction of a 
threshold is deeply connected to the definition of events, 
as indicated in Fig. ^ 

The time series have been downloaded from the 
"Space Physics Interactive Data Resource" (SPIDR) web 
site where each bin represents X-ray flux averaged 
over a given time unit. Among the available signals from 
various GOES satellites H3, we consider the time series 
of the average soft X-ray flux measured in W/m 2 with 
photons in the range from 1 to 8 A. See Table[l]for details. 
We have also isolated two periods, roughly correspond- 
ing to the most recent minimum ( "min" ) and maximum 
("max") of the solar cycle. For comparison, we also an- 
alyze data in publicly available flare catalogs which 
identify events with time intervals when the intensity is 



FIG. 1: (Color online) Explanation of various times used in 
this work. Here the threshold intensity is / = 2 x 10~ 6 W/m 2 . 
According to the standard definition of flares, ai and a? would 
be two separate events. In our case they are separated if, for 
instance, I = 3 x 10~ 6 W/m 2 . This shows that the set of 
events defined by different thresholds are not trivially related 
to each other or to the flares listed in the standard catalogs. 



higher than a local average of the signal. 

In contrast, we define an event, or flare, to be simply 
the interval during which the intensity exceeds a certain 
fixed threshold. To compensate for fixing the threshold 
we study, in detail, the dependence on the threshold value 
and obtain results that are independent of the threshold 
over a wide range. Various times related to our defini- 
tion of events (durations, waiting and quiet times) are 
explained in Fig. ^ One could expect differences when 
events defined by our simple procedure or by the stan- 
dard flare catalogs are analyzed. However, we find that 
the statistical distributions are mostly similar. Hence, 
we expect similar results if other criteria are used to de- 
fine interoccurence times, such as the time difference be- 
tween subsequent maxima in the signal, referred to as 
the "laminar times" in Ref. |8J. Indeed, that definition 
also allows a systematic variation of intensity threshold 
used to select maxima that could be compared with the 
results shown here. 

The number of events with intensity greater than or 
equal to a given threshold I, N(>I), is shown in Fig.^a) 
for the entire data set, at "min" and at "max". For 
the whole catalog the number of flares N(>I) behaves 
approximately as I' 13 , with j3 = 1-2(1) for intensities 
greater than « 1CP 6 W/m 2 . Scaling breaks down be- 
low I w 10 -6 , where N(>I) increases with /. This 
clearly shows that N(>I) is not the cumulative version 
of any probability distribution, because flares at different 
thresholds are different objects. In the two sub-regimes 
"min" and "max", we find power law behavior N(>I) ~ 
/-Anax^ w ith /3 max = 1.2(1), for sufficiently large I. Dur- 
ing the minimum of the cycle, however, another scal- 
ing regime appears. Indeed, N min (>I) ~ 7 _/3min , with 
/5 min = 0.7(1), for I < 1CT 6 W/m 2 . Within statistical 
error, the exponent /? max agrees with the (cumulative) 
distribution of peak fluxes measured by Aschwanden et 
al 0, who obtained /? max + 1 = 2.08 ± 0.03. 
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FIG. 2: (Color online) (a) Number of events with intensity 
greater than a given threshold, for the entire record (*), at the 
minimum (+) and maximum (x) of the cycle. The straight 
lines represent power laws with the quoted exponent, (b) 
The distribution of flare durations, for different thresholds, at 
"min" and "max" of the solar cycle (the former shifted down 
by three units on the log-scale). The symbols are explained 
m Tabled while the curves are fits using Eq. Q. 

The thresholds and associated symbols used by us to 
define the events are given in Table |U For the entire data 
record, we choose five thresholds with / > 10~ 6 , where 
N(>I) is a decreasing function of /. The same thresholds 
are also used in the "max" regime. Since the flux at 
the maximum of solar activity is typically two orders of 
magnitude greater than at the minimum, a definition of 
flares by means of the same set of thresholds could be 
unfeasible. Five different thresholds are used at solar 
minimum to obtain reasonably good statistics. 

We first discuss the distribution of duration times, 
P(td)- It has a power law tail, with critical exponent 
7dur, which extends to longer durations on lowering the 
threshold, as shown in Fig.[2Ib). Previous measurements 



TABLE I: Time series and selected intensity thresholds. 
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"By concatenating signals of satellites GOES 6, 7, 8, and 10 Il7l . 
we reconstruct the time series representing almost two solar cycles. 
Since some data are missing, values in the empty bins are set by the 
last recorded value before each of them. In this way, flat plateaus 
of intensity are created, introducing a minimal bias into the data. 

''Data from GOES-8, minimum of the solar cycle. 

c Data from GOES-10, maximum of the solar cycle. 

^Symbol used in the figures to denote the corresponding data. 
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FIG. 3: (Color online) The distributions of recurrence times, 
progressively shifted along the t/-axes with the threshold, for 
clarity, (a), (c) and (d): Distributions of waiting times, for the 
entire record and at the minimum and maximum of the solar 
cycle, respectively, (b): Distribution of quiet times for the 
entire record. The power law tails are indicated with straight 
lines whose slope is 7. The symbols are explained in Table 
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FIG. 4: (Color online) The rescaled distributions of waiting 
times, arbitrarily shifted vertically to be distinguished. From 
below, for the maximum of the solar cycle, for the minimum, 
for the entire catalog, and for the GOES flare catalog. Fits 
as described in the text are shown for the different phases. 



using a different definition of flares than that put forward 
here found 7d U r = 2.17 to 7d ur = 2.54, depending on 
the range of times fitted , while Litvinenko obtained 
7dur = 2 using dimensional analysis |2(ij|. 

However, P(t d ) crosses over from a power law at large 
times to a constant regime at short times. The entire 
distribution for all thresholds and all time periods is con- 
sistent with the function 

P(t d )~(l +td /t* d )-^ (1) 

with 7d ur = 2.0(1) |2lJ, and t* d ~ 10 min at solar mini- 
mum while 7dur = 2.3(4) and t* d « 20 min at solar maxi- 
mum 22]. 

The waiting and quiet time distributions are shown in 
Fig.[3]for different thresholds and regimes. Each P(t w ) is 
similar to its respective P{t q ), especially at large times, 
where both decay as power laws ~ i~ 7 . Hence, the scale- 
free duration of flares is not giving peculiarities in the 
passage from t q statistics to t w ones. The waiting/quiet 
time exponents 7 at lower thresholds have been evaluated 
in the three regimes. Within statistical error these val- 
ues, indicated in Fig. |3J agree with the ones determined 
by Wheatland and Litvinenko ^(J; wn0 analyzed flare 
catalogs. However, one can observe that by increasing 
the threshold, both P(t w ) and P(t q ) evolve continuously, 
becoming flatter up to longer times for higher thresholds. 
This aspect was not caught by any previous studies [lfj , 
whose results were obtained without systematically vary- 
ing any threshold. 

A scaling argument similar to one recently put forward 
by Bak et al for waiting time statistics of earthquakes [2^ 
can unify in a single scaling function the waiting/quiet 
time statistics. We argue that N(>I) provides the right 



rescaling factor for the recurrence times, namely, the one 
that gives a collapse onto a single curve of all the dis- 
tributions measured with different intensity thresholds, 
/. In particular, we rescale the interoccurrence times by 
their average, which is inversely proportional to the av- 
erage rate of events, R(I) = N(>I)/AT 24], where AT 
is the time span of the record. Thus, the distribution of 
waiting times for a given threshold is given by 

P{t w ,I)~R(I)g(t w R(I)) . (2) 

Unlike the universal waiting time distributions for earth- 
quakes |U 124) . but similar to P(t d ), the scaling func- 
tion g for the flare waiting and quiet times is also well- 
described by the function, g(x) ~ (1 + x/x*)~ J , as shown 
in Fig. 01 Furthermore, the rescaled distributions have 
parameters that depend on the phase of the solar cycle: 
for the whole catalog, 7 ~ 2.16(5) and x* w 0.26. At the 
minimum of the cycle 7 = 1.51(5) (x* w 0.02), while at 
the maximum 7 = 2.83(10) (a;* « 0.85). 

To compare with the standard definition of flares, we 
repeat the above analysis using the GOES flare cata- 
log from 1975-2003 0. Now the threshold / represents 
the peak intensity associated with the flare. The rate of 
events with threshold /, R(I), is the number of events in 
the catalog with peak intensity greater than / divided by 
the total duration of the catalog. Using the same thresh- 
olds as before for the whole catalog, and rescaling the 
distributions using Eq. (01 we obtain the results shown 
in Fig. 0] (upper group of data). In this case power law 
behavior with an exponent 7 = 2.19(5) (x* as 0.28) is 
observed at late times, although at short times there are 
deviations from data collapse at lower thresholds. The 
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deviations may be due to undercounting short waiting 
times following a large flare, an obscuration effect previ- 
ously pointed out by Wheatland [25| . All of the data sets 
can also be fitted with Levy functions, which turn out to 
give comparable results except for the waiting times at 
solar maximum, where the fit with Levy distributions is 
inferior [2(|. We choose here to focus on the fit using a 
single function (Eq. f ) since it appears to describe all the 
data sets equally well. 

Wheatland pjll(J,|25j has modeled the solar flare wait- 
ing time distribution in terms of a time-dependent Pois- 
son process with a flaring rate X(t) . When the flaring rate 
varies slowly over a waiting time, t Wl the distribution of 
waiting times can be written as 



P(t w ,I) 
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F I {X)exp(-Xt w )X 2 dX 



(3) 



where the average flaring rate A/ = J °° Fi(X)\d\ and 
Fj(X)AX is the fraction of the time the rate to produce 
flares exceeding intensity / is within AA of A. The func- 
tion we find to fit the data has //(A) = X 2 Fj(X)/X corre- 
sponding to the Gamma distribution: 
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/ Xx* 
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(4) 



A mathematical equivalence with superstatistics for- 
mulas of Beck and Cohen can be made by mapping 
A — > P, t w — > E and /(A) — > /(/3). No necessary physical 
connection is implied in this equivalence, since one could 
just as well think of a subordination mechanism |2^| as 
being at the basis of Eq. ©. It is also worth remark- 
ing that several turbulent systems have been analyzed 
recently using the superstatistics framework. These in- 
clude velocity differences in Taylor-Couette flow , in- 
termittency of the wind [3(| or solar wind |3l| . 

At solar maximum the critical exponent 7 for the wait- 
ing time distribution is close to 7 = 3. This implies that 
distribution of flaring rates Fj(X) is close to exponential 
for a range of J. Since high intensity flares predominately 
arise from active regions during solar maximum, the ori- 
gin of this distribution could be investigated by tracking 
the flaring rates of individual active regions. 

At solar minimum, the critical exponent 7 « 3/2 for 
a range of intensities /, implying that the distribution of 
flare rates 



FT" (A) ~ A 



-3/2 



exp 



Xx* 
W) 



(5) 



This formula describes the probability distribution for 
the number of offspring in a subcritical branching pro- 
cess |3^ . 

Alternatively, the marginal behavior in on-off intcrmit- 
tency also gives a distribution of laminar times with an 
exponent 7 = 3/2 |l5| . in very good agreement with the 



quiet time distribution at solar minimum. In fact an in- 
termittent on-off dynamo 0] has been used to describe 
the solar cycle and long term records of solar activity 
such as grand minima. Our results lead us to specu- 
late that such a dynamo operating in a marginal state 
may also be able to capture the quiet times of flares - 
excluding active regions. Self-organized criticality may 
provide a mechanism for this dynamo to sustain itself in 
a marginal state. Active regions, superimposed on this 
fluctuating state, may represent plasma instabilities in 
the on-off dynamo, with their own emergent behavior. 
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